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We apply the Schrodinger Functional (SF) formalism to determine the renormalisation group running of four- 
fermion operators which appear in the effective weak Hamiltonian of the Standard Model. Our calculations are 
done using Wilson fermions and the parity-odd components of the operators. Preliminary results are presented 
for the operator Ova = (s7^d) (s7^75d). 



1. Introduction 

Weak matrix elements (WMEs) of parity-odd 
four-fermion operators are crucial in many Stan- 
dard Model processes (e.g. A/ =1/2 enhance- 
ment, AS* = 1 non-leptonic Kaon decays). The 
non-perturbative aspects of these decays are ex- 
pressed as WMEs {Tnr\Ok\K) of parity-odd oper- 
ators Ok (where, in standard notation, k — VA± 
AV, SPztPS, TT denotes the various Dirac struc- 
tures). On the other hand, the non-perturbative 
aspects of K° — K° oscillations (i.e. AS* — 2 tran- 
sitions) are studied through WMEs of parity-even 
operators {K°\Ovv+aa\K^) ■ Even in this case, 
however, it is possible to map the problem into 
the computation of parity-odd matrix elements, 
either through the use of Ward Identities 0] or 
by implementing the twisted mass QCD formal- 
ism [||. 

Unlike their parity-even counterparts, all 
parity-odd operators of interest display renor- 
malisation patterns which are unaffected by the 
breaking of chiral symmetry on the lattice (i.e. 
the Wilson term) [0. Our program consists in 



computing non-perturbatively the anomalous di- 
mension (AD) matrix of the complete basis of 
parity-odd dimension-six operators with Wilson 
fermions. Our renormalisation scheme of choice 
is the SF, which is mass independent. This has 
two implications. First, our results also apply 
to those operators which characterise A/ =1/2 
transitions. Although the complete removal of di- 
vergences in this case requires complicated power 
subtractions, these, being proportional to quark 
masses, do not affect the operator AD. Sec- 
ond, our results also apply to operators of heavy 
flavour transitions (e.g AC, AS = 2, AB = 1). 

The renormalisation of any four-fermion oper- 
ator (up to power subtractions) can be mapped 
to that of an operator with the same colour-Dirac 
structure but four distinct flavours. For example, 
once the operators 
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+ (V'i7m75'(/'2) (V'37/.V'4) ± (2 ^ 4)] (1) 

have been properly renormalised, we can map the 
four distinct flavours to the physical flavours con- 
cerning the process of interest. Here we will only 
present results for the operator Q = Oyj^^^y. 
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2. The step scaling function (SSF) 

Our work consists in a straightforward gener- 
alisation of the computation of the SSF of the 
quark mass, performed by the Alpha collabora- 
tion in Q . The operator anomalous dimension 7 
is defined through 



dQ 



7(5) Q 



{Q and g are renornialised quantities), 
turbative expansion is given by 



(2) 
Its per- 



7(5) -9^ (70 + 71.9^ + 725"* + ■ • ■) , (3) 

with 7o — l/(47r^) a universal coefficient. The 
(scheme-indepedent) renormalisation group in- 
variant (RGI) operator is obtained upon integrat- 
ing eq. (|) 
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(60 is the LO coefficient of the Callan-Symanzik 
/3-function) . Once the scale dependence of the 
renormalised coupling g in the SF scheme is 
known , and proceeding along the line of ref . , 
three elements are required for the complete NP 
determination of a WME of the form {f\Q^^^\i)'- 
(i) the lattice bare matrix elements at several 
bare couplings (work in progress); (ii) the value 
of the renormalisation constant Zq at the same 
bare couplings and a fixed low-energy reference 
scale of the order of Aqcd', (iii) the operator SSF 
from this non-perturbative low-energy scale to 
high ones, where g{fi) is safely in the perturbative 
region and the integral in eq. can be performed 
analytically to determine the ratio Q^'^^ /Q{lJ.)- 

We regularise the theory on a lattice of physical 
size with standard SF boundary conditions . 
Renormalisation is performed in the chiral limit, 
at a scale equal to the IR cutoff L~^. Since the 
operator renormalises multiplicatively, 

Q(i-i)= hm Zgigo^L/a) Q{a) , (5) 

a — >u 

the UV divergence is removed by a single renor- 
malisation condition, namely by setting the ex- 
pression 

/i~^^^(5o,i/a) ZQ{ga,L/a) hgixo) (6) 



at xq = L/2 to its tree level value. The corre- 
lator hg is as shown in Fig. |l]; /i serves to can- 
cel the quark boundary field renormalisation (see 
ref. 0). For more details and an explicit nota- 
tion, see ref. |6). Different choices of boundary 




Figure 1. The SF correlation function hq. 
Note that the "spectator" quark line connecting 
the two time boundaries corresponds to a fifth 
flavour. 



operators correspond to different renormalisation 
schemes. Parity conservation and the properties 
of SF boundary quark fields limit us to five non- 
trivial independent choices, 

75 ^ 75 75 75 ^ 7fe 7fe 
7fc ^ Ik 75 Ik ^ 75 7fc 

li ^ 1] Ik <^ijk , (7) 

plus their time-reversed counterparts (in eq. (^ 
the 7-matrix on the Ihs. of the symbol <-> cor- 
responds to the Dirac structure at time xq — 0, 
whereas the two 7-matrices on its rhs. correspond 
to the Dirac structures at time = L). 
The SSF is defined as 

CTQ u) = Imi Eq M,a/L = 1™ TTT ^ 

a^o a^a ZQ(go,L/a) 

(u = g^{L)) and, taking into account the run- 
ning coupling from ||^, can be used to run renor- 
malised matrix elements along a broad range of 
scales through the recursion 



(9) 



Our quenched simulations have been perfomed at 
the parameters (volumes, couplings and hopping 
parameters) of ref. [Q . The most difficult numer- 
ical task consists in the extrapolation of our data 
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Figure 2. Continuum limit extrapolation of 
the SSF at renormalised couplings u — 1.2430 
(upper plot) and u = 2.770 (lower plot). Solid 
(open) points are data corresponding to simula- 
tions with an improved Clover (Wilson) action. 
The leftmost point in each plot is the result of a 
combined linear extrapolation to the continuum 
limit of both data sets, with slopes indicated by 
the two lines. In both cases, the data at largest 
lattice spacing have not been included in the fit. 



to the continuum limit: although our action is 
the non-perturbatively improved Clover one (i.e. 
cutoff effects are O{o?) from the bulk, and ©(^qo) 
from the time boundaries) , our operator is unim- 
proved, which gives rise to 0[a) effects. Thus, for 
some renormalised couplings the continuum limit 
extrapolation appears to be stable, whereas for 
others it is rather unreliable. In order to improve 
on this we are currently repeating our computa- 
tion, in the spirit of ref. Q with an unimproved 
Wilson action. Since the SSF corresponding to 
the same renormalised coupling has the same con- 
tinuum limit, combined fits of the Clover and Wil- 
son SSF should stabilise the results. To show that 
this is actually the case, we present in Fig. || two 
of our first combined extrapolations at two dif- 
ferent values of the renormalised coupling. Until 
this analysis is completed, we shall refrain from 
quoting an estimate for the quantity Q^^^ /Q{fJ.)- 
In Fig. ^ we show our present results with an 
improved action for the SSF associated with the 
first renormalisation condition in eq. (M). The 
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Figure 3. Step scaling function ctq in the contin- 
uum limit. The solid line is the LO perturbative 
prediction for erg. 



step scaling functions associated with the other 
choices of renormalisation scheme are qualita- 
tively similar. Our results need to be supple- 
mented by the NLO calculation of the SSF, which 
will enable us to estimate the importance of NP 
contributions to this quantity in the SF scheme. 
It will also provide the necessary information for 
the matching between SF and MS schemes. This 
calculation is well under way. 
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